Enhanced photon production from quark-gluon plasma: Finite-lifetime effect 
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Photon production from a thermalized quark-gluon plasma of finite lifetime is studied directly in 
real time with a nonequilibrium formulation that includes off-shell (energy nonconserving) effects. To 
lowest order we find that production of direct photons form a quark-gluon plasma of temperature T ~ 
200 MeV and lifetime t ~ 10 — 20 fm/c is strongly enhanced by off-shell (anti)quark bremsstrahlung 
q(q) — > q(q)^y. The yield from this nonequilibrium finite-lifetime effect dominates over those obtained 
from higher order equilibrium rate calculations in the range of energy E > 2 GeV and falls off with 
a power law for E 3> T. 
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The observation of a novel phase of matter, the quark- 
gluon plasma (QGP), is one of the most important goals 
of ultrarelativistic heavy ion experiments currently un- 
dertaken at CERN SPS and BNL RHIC @]. The quark- 
gluon plasma formed in the early stage of the collision ex- 
pands and cools rapidly to a mixed phase of quarks, glu- 
ons, and hadrons, and ultimately undergoes a freeze-out 
from a state of hadronic gas. Estimates based on energy 
deposited in the central collision region at RHIC energies 
•y/s ~ 200^4 GeV suggest that the lifetime of a deconfincd 
phase of quark-gluon plasma is of order 10 — 20 fm/c with 
an overall freeze-out time of order 100 fm/c g|. An im- 
portant aspect is an assessment of nonequilibrium effects 
associated with the rapid expansion and finite lifetime 
of the plasma and their impact on experimental observ- 
ables. 

Amongst various experimental signatures proposed to 
detect the quark-gluon plasma phase, photons (both 
real and virtual) have longbeen considered as the most 
promising direct signals J5|. This is because, unlike 
strongly interacting hadrons, photons have a mean free 
path much larger than the typical size of the plasma 
formed in ultrarelativistic heavy ion collisions. Once pro- 
duced they escape from the system without further in- 
teraction, thus carrying clean information from the early 
hot quark-gluon plasma phase. 

The first observation of direct photon production in 
ultrarelativistic heavy ion collisions has been recently re- 
ported by the CERN WA98 Collaboration in Pb+Pb col- 
lisions at sfs = 158 A GeV Q|. The transverse momen- 
tum distribution of direct photons is determined on a 
statistical basis and compared to the background photon 
yield predicted from a calculation of the radiative decays 
of hadrons. The most interesting result is that a signifi- 
cant excess of direct photons beyond that expected from 
proton-induced reaction at the same y/s is observed in 



the range of transverse momentum greater than about 
1.5 GeV/c in central collisions. While it is not yet clear 
whether a QGP was formed in the central collision region 
at SPS energies, this result does suggest the experimental 
feasibility of direct photon production as a signal of the 
QGP phase, expected to be formed at RHIC energies. 

One goal of this article is to study directly in real 
time the effect of the finite QGP lifetime on direct pho- 
ton production. We focus on the direct photon yield 
from a thermalized quark-gluon plasma of temperature 
T ~ 200 MeV and lifetime t ~ 10 - 20 fm/c in accor- 
dance with estimates based on collision energies reached 
at RHIC. Another goal is to compare this nonequilib- 
rium photon yield to those of previous investigations (g|Jfj] 
that obtain an equilibrium production rate by taking the 
thermal average of transition amplitudes, which in ef- 
fect is tantamount to assuming an infinite lifetime for 
the QGP. It is not the purpose of this article to compare 
direct photon production from the QGP with that from 
the hadronic gas, as this has been studied in detail in 
Refs. if). 

Many investigations in the literature have been de- 
voted to hard real photon production from the quark- 
gluon plasma j5|-|8|. Assuming a QGP in thermal equi- 
librium with an infinite lifetime, these authors focused on 
the equilibrium photon productionj^atCj which for pho- 
tons of momentum p is given by 



E 



dN 



d 3 pd 4 x (2tt) : 



■n B {E)lmIL R (E). 



(1) 



Here E = |p|, H R (E) is the retarded transverse photon 
self-energy at finite temperature T evaluated on shell, 
and ns{E) = \/{e E / T — 1) is the Bose-Einstein distri- 
bution. Using the hard thermal loop (HTL) resummed 
effective perturbation theory developed by Braaten and 
Pisarski |L0| , Kapusta et al. || and Baier et al. || showed 
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that at one-loop order (in effective perturbation theory) 
the processes that contribute to photon production are 
the gluon-to-photon Compton scattering off (anti)quark 
9(9)5 ~~ * 9(9)7 and quark- antiquark annihilation to pho- 
ton and gluon 99 — ► g-j. The corresponding rate of ener- 
getic (E 3> T) photon emission for two light quark flavors 
(u and d quarks) is given by (5| 
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d 3 pd 4 x 



one-loop 
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where a is the fine-structure constant and a s — gl/^K 
with g s being the strong coupling constant. In a recent 
development Aurenche et al. [f| have found that the two- 
loop contributions to the photon production rate arising 
from (anti)quark bremsstrahlung 99(3) — > 99(3)7 and 
quark-antiquark annihilation with scattering qqq(g) — > 
9(3)7 are °f the same order as those evaluated at one 
loop. The two-loop contributions to the photon produc- 
tion rate read H 
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d 3 p d A x 



two-loop 
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where Jt ~ 4.45 and Jl ~ —4.26 for two light quark fla- 
vors. Most importantly, they showed that the two- loop 
contributions completely dominate the photon emission 
rate at high photon energies M. 

To study photon production from a QGP of finite 
lifetime, we use a real-time kinetic approach ||il| based 
on noncquilibrium quantum field theory |I2j-|l4||, which 
when improved by a resummation via a dynamical 
renormalization group provides a consistent microscopic 
derivation of quantum kinetics from the underlying the- 
ories |i||L4|] • One of the advantages of this real-time ki- 
netic approach is that it is capable of capturing off-shell 
(energy nonconserving) effects originating in the finite 
system lifetime, as completed collisions are not assumed 
a priori (l4| . 

Because of the abelian nature of the electromagnetic 
interaction, we will work in a gauge invariant formulation 
in which physical observables (in the electromagnetic sec- 
tor) are manifestly gauge invariant pIR ■ Since photons 
escape directly from the quark-gluon plasma without fur- 
ther interaction, it is adequate to treat them as asymp- 
totic particles. In the Heisenberg picture the number 
operator N(p, t) that counts the total number of photons 
of momentum p (per phase space volume) at time t is 
defined by 



N(p,f)=^4(p,t)a A (p,<), 



(4) 



A=l 



where a\(p, t) [a\{p, t)] is the annihilation (creation) op- 
erator that destroys (creates) a photon of momentum p 



and polarization A at time t. The time- dependent pho- 
ton production rate, which up to a trivial factor E/(2ir) 3 
is related to the expectation value of the time deriva- 
tive of N(p,t), can be obtained by using the Heisenberg 
equations of motion (for details, see Ref. [3). In the 
framework of nonequilibrium quantum field theory, the 
time-dependent photon production rate is given by |14[] 
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(5) 

where k = p + q. Here e q is the electromagnetic cou- 
pling constants of the quarks, and tp denotes the (gauge 
invariant) quark field. The "+" ("— ") superscripts for 
the fields refer to fields defined in the forward (backward) 
time branch |l2[ . 

As usual we assume that there are no photons initially 
and those once produced will escape from the plasma 
without building up their population. Therefore the 
QGP in effect is treated as the vacuum for the photons. 
Consequently, the nonequilibrium expectation values on 
the right-hand side of Eq. (g|) is computed perturbatively 
to order a and in principle to all orders in a s by using 
real-time Feynman rules and propagators. Furthermore, 
the photon propagators are the same as those in the vac- 
uum. 

We shall further assume the weak coupling limit a <C 
a s <C f • Whereas the first limit is justified and essen- 
tial to the interpretation of electromagnetic signatures 
as clean probes of the QGP, the second limit can only be 
justified for very high temperatures, and its validity in 
the regime of interest can only be assumed so as to lead 
to a controlled pcrturbative expansion. 

It is not hard to see that the lowest order contribu- 
tion to Eq. (J5J) is of one quark loop and of order a. In 
using bare quark propagators we consider the quark mo- 
mentum in the loop to be hard, i.e., 9 > T. Soft quark 
lines require HTL resummed effective quark propagator 
leading to higher order corrections. Indeed, the one-loop 
diagram with soft quark loop momentum is part of the 
higher order contribution of order aa s that has been cal- 
culated in Refs. [||||. 

In this article we focus on the lowest order contribu- 
tion which, as will be understood below, contributes to 
direct photon production solely as a consequence of the 
finite QGP lifetime, and has therefore been missed by all 
previous investigations which assumed an infinite QGP 
lifetime. 

For two light flavors (u and d quarks), the lowest or- 
der time-dependent photon production rate is found to 
be given by H] 



E 



dN{t) _ 2 
d 3 pd 4 x ~ (2tt) 3 



dw1Z(bj) 



9in[(u> - E)(t - to)] 



it(uj - E) 



(6) 



with 



n(w) = 



20n 2 a f d 3 q r r , A - w , AN 

x n F (fc)[l — 7ii?(g)]<5(o; — k + q) 
+ [l + (p-k)(p-q)]{n F (fc)n F (g) 
xj(w-ife-g) + [l- n F (fc)] 
x[l-n F ( g )]5(w + fc + g )}}, 



(7) 



where t is the initial time at which the QGP is formed, 
q = |q|, and np(q) = \j{e q l T + 1) is the Fermi-Dirac 
distribution. 

A detailed analysis shows that the first delta function 
S(u> — k + q) with support below the light cone (u> 2 < E 2 ) 
corresponds to the Landau damping cut, and the last 
two delta functions 5(lu =f k =F q) with supports above 
the light cone correspond to the usual two-particle cut. 
Furthermore, one recognizes that 1Z(lu) has a physical 
interpretation in terms of the following off-shell (energy 
nonconserving) photon production processes: the first 
term describes (anti)quark bremsstrahlung q(q) —> q{q)^] 
the second term describes quark-antiquark annihilation 
to photon qq — > 7; and the third term describes creation 
of a photon and a quark-antiquark pair out of the vac- 
uum — > qq'y. We would like to stress that the photon 
and (anti)quark in these processes are all on-shell parti- 
cles, and off-shellness refers to nonconservation of energy 
as will be explained below. It is also noted that the first 
two processes are inherently distinct from the third one. 
Whereas the first two processes arise solely due to the 
presence of the QGP, the third one in effect is a "vac- 
uum" process that is Pauli suppressed during the QGP 
lifetime. 

It can be shown that lZ(u>) is related to the imagi- 
nary part of the one- loop retarded transverse photon self- 
energy with hard loop momentum p3 



TZ(uj) = -ns(w)ImII i? (w)| onc _ lo0 p. 



(8) 



If, as is usually done, the QGP is assumed to be of in- 
finite lifetime, then we can take the infinite time limit 
to — > —00 in the argument of the sine function in Eq. (1q) 
and use the approximation 



sin[(w - E)t] t- 

tt(w - E) 



5(u> - E). 



(9) 



This is the assumption of completed collisions that is in- 
voked in time-dependent perturbation theory leading to 
Fermi's golden rule and energy conservation. Under this 
assumption Eq. (fy) is recovered to lowest order and one 
finds a time-independent photon production rate propor- 
tional to 1Z(E), provided that the latter is finite. 

In the present situation, however, the three delta func- 
tions in lZ(u>) cannot be satisfied on the photon mass 
shell. Therefore under the assumption of completed col- 
lisions the off-shell contribution to the photon produc- 
tion rate simply vanishes due to kinematics. Physically 



this reflects the energy nonconserving nature of the cor- 
responding photon production processes. The "vacuum" 
process is independent of the presence of the QGP or its 
lifetime and persists for a infinitely long time. There- 
fore for the third term of TZ(u>) we have to take the infi- 
nite time limit, which leads to the vanishing of the off- 
shell "vacuum" process — > qq^f by energy conservation. 
Only the "medium" processes depend on the presence 
and finite lifetime of the QGP, hence we only need to 
consider off-shell (anti)quark bremsstrahlung and quark- 
antiquark annihilation in the rest of our discussion. 

For any finite QGP lifetime the time-dependent rate 
given in Eq. (ra) is finite and nonvanishing, thus leading 
to a nontrivial contribution to direct photon production. 
In this article we focus on the photon yield instead of the 
photon production rate, since the former is the quantity 
of phenomenological interest that can be measured in ex- 
periments. The photon yield is obtained by integrating 
the rate over the lifetime of the QGP. Using Eq. (p), one 
obtains 



E 



dN(t) _ 2 
d 3 pd 3 x ~ (2tt) 3 



dujTZ(uj) 



1 - cos[(w - E)t] 



tt(lu - E) 2 



(10) 



where, here and henceforth, we have set to = 0. 

Before proceeding to a numerical study, we give an 
analytic estimate of the behavior of the photon yield in 
the HTL approximation. In this approximation the lead- 
ing contribution to Imll (o>) | on c-ioop is dominated by the 
Landau damping cut P,H, which corresponds to off- 
shell (anti)quark bremsstrahlung. Thus, from Eq. (jg), 
we find || 



"KhtlM 



20 tt 2 aT 2 uj 
T 12 E 



to 
E 2 



n B (uj)9(E 2 -uj 2 ). 
(11) 



(12) 



For £<T, 7£htl(w) can be further simplified as 

m , 1 b«t20 7t 2 «T 3 / u 2 \ n ,„ 2N 



The dominant contribution of the w-integral in Eq. (Il 
for E <C T arises from the region where the resonant de- 
nominator vanishes, i.e., u w E [13 14 1. Using Eq. (|l 
we obtain |14] 



dN(t) s<t 5 a T 3 , 



(13) 



where 7 = 0.577... is the Euler-Mascheroni constant. 
It has been shown in Ref. Q that for photons of en- 
ergy E ^ T the finite-lifetime contribution to the photon 
yield, Eq. dlq), is comparable at early times to that of 



order aa s obtained from the equilibrium rate given in 
Eq. © 

Fig. [l] shows that for E ^$> T, TZwtl{^>) is exponen- 
tially suppressed in the region T < lu < E. From this 
observation we emphasize that (i) because for E 3> T the 
threshold contribution near the photon mass shell lu k E 
is exponentially suppressed, the w-integral in Eq. (fill) is 
now dominated by the interval — E < lu < T, which cor- 
responds to highly off-shell (anti)quark bremsstrahlung; 
(ii) as the integrand in Eq. ( |l0|) is positive-definite and 
the function f — cos[(cj — E)t] averages to 1 for large 
t in the region —E < lu < T, for fixed E > T the 
yield approaches a constant at large times; and (iii) in 
contrast to those obtained from equilibrium rates, the 
yield for E 3> T is not suppressed by the Boltzmann 
factor e~ E ' T . These important nonequilibrium aspects 
have noteworthy phenomenological implications studied 
in detail below. 

We now perform a numerical analysis of the nonequi- 
librium finite-lifetime contribution to photon produc- 
tion directly in terms of the yield given by Eqs. ( |I(]| ) 
and (in) [but, as explained above, without including the 
third term of 1Z{lu)] and compare the results to those 
obtained from the equilibrium rates given in Eqs. (0) 
and (||). For this we consider a QGP of temperature 
T = 200 MeV and lifetime t = 10 - 20 fm/c. As re- 
marked above, this is approximately the scales of the 
QGP temperature and lifetime expected at RHIC ener- 
gies. For the value of a s at this temperature, we use 
a s (T) = 6tt/(33 - 2JV» ln(8T/T c ) ||, where N f = 2 is 
the number of quark flavors and T c is the quark-hadron 
transition temperature. Taking T c = 160 MeV from 
lattice QCD calculation for two quark flavors, we find 
a s w 0.3 at T = 200 MeV. 

The results of the numerical study of the nonequilib- 
rium photon yield are displayed in Figs, g and ||. Fig. g 
depicts the time evolution of the nonequilibrium yield up 
to 20 fm/c for photons of energies E = 1 and 2 GeV. 
We observe that the yield exhibits clearly a logarith- 
mic time-dependence at late times during the QGP life- 
time. Numerical evidence shows that the nonequilibrium 
yield EdN(t)/d 3 pd 3 x falls off with a power law E~ v with 
v k, 2.14 for E ^> T. In addition, the numerical result 
also reveals that the dominant photon production process 
is (anti)quark bremsstrahlung q(q) — > q(q)"f. 

Figured shows a comparison of the nonequilibrium and 
equilibrium (one-loop and two-loop) contributions to the 
direct photon yield in the range of energy T < E < 
4 GeV. Whereas the two-loop contribution dominates 
the direct photon yield for smaller values of E, a sig- 
nificant enhancement of the direct photon yield due to 
the nonequilibrium contribution is seen at E > 2 GeV 
as a consequence of its power law falloff for E ^> T. 
In particular, the nonequilibrium contribution is larger 
than the equilibrium contributions by several orders of 
magnitude for E > 3 GeV. As the nonequilibrium con- 
tribution for fixed E 3> T approaches a constant at large 
times, the equilibrium contributions, which grow linearly 



in time, will eventually dominate the yield if the QGP 
has a very long lifetime. However, the linear growth in 
time of the equilibrium contributions has to compensate 
the Boltzmann suppression for E 3> T. Therefore we 
emphasize that for E 3> T the equilibrium contributions 
could dominate the yield if the QGP lifetime is of or- 
der 10 2 — 10 3 fm/c or larger, which nevertheless is very 
unrealistic at RHIC energies. 

Our analysis can be extended to the case of direct pho- 
ton production from a QGP away from chemical equilib- 
rium by replacing the equilibrium (anti)quark distribu- 
tions with the undersaturated ones [|6[. We expect the 
main features to remain in the more general situation. 

In this article we have studied directly in real time the 
production of direct photons from a thermalized QGP of 
temperature T ~ 200 MeV and lifetime t ~ 10-20 fm/c. 
To lowest order in perturbation theory we find that di- 
rect photon production features a power law spectrum 
for E ^S> T due to a significant enhancement by off-shell 
(anti)quark bremsstrahlung. 

Our main conclusion is that the novel nonequilibrium 
effect originating in the finite lifetime of a thermalized 
QGP leads to an enhancement in the direct photon yield 
for E > 2 GeV, which dominates over the yields obtained 
from equilibrium rate calculations for the same QGP life- 
time. To establish a direct contact with the experimental 
observations, the next step of our program will be to in- 
clude effects of hydrodynamical expansion of the QGP 
as well as an assessment of photon production from a 
hadronic gas with a finite lifetime, from hadronization to 
freeze-out. Work in this direction is in progress. 
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FIG. 2. The nonequilibrium yield from a QGP of temper- 
ature T = 200 MeV is plotted as a function of time for E = 1 
(upper) and 2 GeV (lower). The circles denote the numerical 
result and the solid line is a logarithmic fit. 
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FIG. 3. Comparison of various contributions to the di- 
rect photon yield from a quark-gluon plasma of temperature 
T = 200 MeV and lifetime t = 10 fm/c. The inset shows the 
figure on a log- log plot. 



FIG. 1. The function 7?.htl(^) is plotted in the limit 
E > T. Here we take E = 2 GeV and T = 200 MeV. 



